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ICASE Fluid Mechanics 


Due to increasing research being conducted at ICASE in the field of fluid mechanics, 
future ICASE reports in this area of research will be printed with a green cover. Applied 
and numerical mathematics reports will have the familiar blue cover, while computer science 
reports will have yellow covers. In all other aspects the reports will remain the same; in 
particular, they will continue to be submitted to the appropriate journals or conferences for 
formal publication. 
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SOME ASPECTS OF THE AEROACOUSTICS 
OF HIGH-SPEED JETS 1 


Sir James Lighthill 
Mathematics Department 
University College London 
Gower Street 
London WC1E 6BT 
ENGLAND 

ABSTRACT 

The Lecture begins by sketching some of the background to contemporary jet aeroacous- 
tics. Then it reviews scaling laws for noise generation by low-Mach-number airflows and 
by turbulence convected at “not so low” Mach number. These laws take into account the 
influence of Doppler effects associated with the convection of aeroacoustic sources. 

Next, a uniformly valid Doppler-effect approximation exhibits the transition, with in- 
creasing Mach number of convection, from compact-source radiation at low Mach numbers 
to a statistical assemblage of conical shock waves radiated by eddies convected at supersonic 
speed. In jets, for example, supersonic eddy convection is typically found for jet exit speeds 
exceeding twice the atmospheric speed of sound. 

The Lecture continues by describing a new dynamical theory of the nonlinear propagation 
of such statistically random assemblages of conical shock waves. It is shown, both by a 
general theoretical analysis and by an illustrative computational study, how their propagation 
is dominated by a characteristic “bunching” process. That process - associated with a 
tendency for shock waves that have already formed unions with other shock waves to acquire 
an increased proneness to form further unions - acts so as to enhance the high-frequency 
part of the spectrum of noise emission from jets at these high exit speeds. 


'This is the original version of the Theodorsen Lecture which was supported by the Institute for Computer 
Applications in Science and Engineering (ICASE) and the National Aeronautics and Space Administration 
(NASA), Langley Research Center, Hampton, VA 23681-0001 under NASA Contract No. NASI- 19480. 
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1. Introduction 


I warmly appreciate the invitation to give the inaugural Theodorsen Lecture in 
honor of a renowned Langley scientist, Theodore Theodorsen. It is moreover a 
very special pleasure for me to pay tribute here to the deep and intricate 
aerodynamic researches of Dr Theodorsen because they powerfully influenced all 
work on airfoil design during the early 1940s and, in particular, some extensive 
activities in this field pursued by the team in the NPL Aerodynamics Division, 
under Sydney Goldstein’s leadership, of which I then formed part. This team’s 
achievements were later to be comprehensively expounded by Dr Goldstein in his 
Wright Brothers Lecture (Goldstein 1948), where he went out of his way to stress 
the fundamental importance of Theodore Theodorsen’s key contributions to airfoil 
theory and to airfoil design. 

By that time, in the late 1940s, I myself was starting work on a rather new 
branch of aerodynamics, which I then called the study of sound generated 
aerodynamically (Lighthill 1952) and which came later to be known as 
aeroacou sties. From the outset a vitally important branch of aeroacoustics was 
concerned with the noise of turbulent jets, a subject to which I was to devote my 
own Wright Brothers Lecture (Lighthill 1963). 

Fundamental investigations of turbulent flows had, of course, been yet 
another of the major interests pursued by Theodore Theodorsen. I have accordingly 
felt that a lecture concerned with turbulent jets and with one of the gravest 
environmental problems posed by their use for aircraft propulsion purposes 
namely, the noise they generate — would form a fitting memorial tribute to that 
great aerodynamicist. 

Furthermore, a lecture on this theme is particularly appropriate to an 
important occasion here at Langley; where early pioneering researches on 
aeroacoustics, including some crucial experiments on jet noise, were carried out 
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by Harvey Hubbard and others (Lassiter & Hubbard 1952). Last year, moreover, 
ICASE and NASA Langley Research Center joined forces to host an exceptionally 
fine Workshop, which gave an authoritative view of the past, present and fiiture 
of aeroacoustics, with special emphasis on Computational Aeroacoustics. This 
Workshop looked forward towards "a second golden age of aeroacoustics", in 
which both new theoretical and new computational approaches, closely linked to 
new experimental techniques, would be employed in response to the challenge of 
those toughened environmental requirements that are now to be imposed on all new 
aircraft designs (Hussaini & Hardin 1993). 

I personally was delighted to participate in this workshop by giving an 
introductory lecture looking back to some of the earlier fundamental discoveries 
in aeroacoustics, as well as by chairing the Final Panel Discussion which looked 
forward to likely future developments as aeroacoustics enters its second golden 
age. In between, I was stimulated by hearing many brilliant contributions 
concerned with meeting some of those new challenges to which I just referred. 

One of the most exciting of these is posed by the US High Speed Civil 
Transport project (HSCT), an extremely promising plan for a supersonic transport 
aircraft ingeniously designed to minimise the level of supersonic-boom annoyance. 
However the corresponding problems of reducing engine noise for such an aircraft 
to within acceptable limits raise some thorny questions and may, in particular, 
demand that a fundamental aeroacoustic analysis of jets at relatively high speed be 
undertaken. 

In a purely aeroacoustical context the appropriate definition of the Mach 
number M of a jet is the ratio of its exit speed U to the speed of sound c a in the 
atmosphere into which it is radiating. Now, the general trend of aeroacoustic Mach 
numbers for civil aero-engines (in other words, for the engines of those aircraft 
which face the greatest aeroacoustic challenges) has been a downward trend for 
very many years, and this has allowed engines to become simultaneously quieter 
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and more powerful because in a wide range of Mach numbers (see figure 1 below) 
the ratio tj of acoustic power radiated to jet power delivered varies as M 5 .At lower 
jet Mach numbers, then, noise radiated can be less even though jet power is 
greater. 

Evidently, it has only been the introduction of aeroengines of ever larger 
and larger diameter that has permitted such increased engine powers to be 
delivered at the necessary low noise levels by means of this downward trend in jet 
Mach number. For a supersonic aircraft, on the other hand, the use of extra-wide 
engines is out of the question because their supersonic-boom emissions, and also 
the associated shock-wave drag, would be unacceptably great. Such considerations 
rule out any similar reduction in jet Mach number in this case. 

It follows that work on engine-noise reduction for the HSCT project has 
been calling for many new fundamental studies of the generation of noise by jets 
at relatively high Mach numbers. This is the range of values of M in excess of 2 
for which (see figure 1 below) the proportion rj of jet power that is converted into 
noise approaches an asymptotic value of 0.01 or a little less. Moreover the noise 
field has become highly directional, because the turbulent eddies that generate 
sound are themselves being convected at supersonic speeds so that they emit their 
own supersonic booms in the Mach direction defined by an eddy convection 
velocity and the atmospheric speed of sound (Ffowcs Williams & Maidanik 1965). 

In an actual aero-engine installation, of course, there may be a great 
difference between the pure noise field of the high-speed jet itself and the overall 
noise radiated from the installation taken as a whole. However, an essential pre- 
requisite for designing that installation so as to bring down ground noise levels is 
to understand as well as possible the primary noise field generated by the high- 
speed jet. This is why so many stimulating lectures contributed to last year’s 
Workshop at Langley (Hussaini & Hardin 1993) were devoted to different aspects 
of that primary noise generation. 
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In this Theodorsen Lecture, the main material which I shall present 
describes certain new researches into which I was drawn as a direct result of 
listening to all those stimulating contributions which, as I have indicated, were 
concerned with how the high-speed jet generates a noise field emitted largely in the 
Mach direction. In fact the qualifying phrase which appears in my title "Some 
aspects of the aeroacoustics of high-speed jets" is mainly intended to recognize that 
noise-generation mechanisms in such jets have long been studied and, above all, 
are being investigated actively today. Accordingly I shall refer only briefly to those 
generation aspects of the problem and concentrate rather on some other aspects 
which appear, relatively speaking, to have been neglected. 

These aspects of the aeroacoustics of high-speed jets on which I shall 
concentrate are related to those effects of nonlinear sound propagation which 
immediately start to modify the noise field once it has been generated. Such 
nonlinear propagation effects may readily be expected to be important from that 
analogy with supersonic-boom generation which I already mentioned; and which 
tends to suggest that each supersonically convected eddy will emit in the Mach 
direction a boom-like signal that should include one or more conical shock waves. 
Accordingly a random sequence of eddies should generate a random, thickly 
packed assemblage of conical shock waves, and this idea is supported by various 
experimental, theoretical and computational studies. 

My lecture, then, is primarily concerned with the nonlinear acoustic 
propagation of such random assemblages of conical shock waves once they have 
been formed. Thus it includes a study of the inherent tendency of the shock 
strengths to decay (through conical spreading and internal dissipation), as opposed 
by increases in strength that occur whenever adjacent shock waves unite. A certain 
"bunching" tends to arise, because a union of two adjacent shocks is found to 
increase the likelihood of further union with other neighbouring shocks. The high- 
frequency part of the noise spectrum is made more intense than would otherwise 
be the case by these bunching tendencies. 
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A completely general theoretical analysis is used to demonstrate the 
universal tendency towards bunching. Then a computational study is carried out in 
order to exemplify details of the process, which is expected to be important for the 
appreciation of how Mach wave fields are modified in the region of conical 
propagation that surrounds a high-speed jet. 

All of this new material concerned with nonlinear propagation is preceded, 
however, by a simplified summary account of how jet noise is generated. This 
account — broadly along the lines of my Wright Brothers Lecture — shows how 
some sort of continuous transition can be discerned between more familiar 
processes of jet-noise generation at relatively low Mach number and that generation 
of a random assemblage of conical shock waves which sets in at the higher jet 
Mach numbers. 

Next, a well known transformation of coordinates (Whitham 1956, Lighthill 
1978) is used to reduce the problem of conical propagation, concerned with how 
a temporal waveform varies with distance r along the Mach direction, into a plane- 
wave problem. In this latter problem the time t replaces r' A , and attention is 
focussed on how a given spatial waveform varies with time. 

The spatial waveform that needs to be studied in this latter context turns out 
to be a spatially unlimited assemblage of random "sawtooth” waves, of the type 
(see figure 2 below) into which a general plane sound wave of large amplitude 
would evolve during a certain time t. It consists of randomly located shocks (with 
random strengths) separated by expansion waves in which the slope of excess 
signal velocity as a function of position takes the value 1 It. 

I shall present a new and quite general theory of the nonlinear dynamics of 
random sawtooth waves. It shows how the inherent tendency of the shock strengths 
themselves is to decay like 1/r but that this tendency is opposed whenever shocks 
unite with adjacent shocks. Moreover those "bunchings" of unions to which 
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I already referred combine in a sort of "snowball" effect to make major 
modifications in the waveform. 

After a transformation back into the original variables describing conical 
propagation, t is replaced by r' h but predicted shock strengths have to be divided 
by a further r h factor. Thus the inherent tendency of the shock strengths is now 
to decay like Mr but we shall yet again find that this tendency is powerfully 
opposed by those "bunching" effects whose special relevance, in a noise context, 
is to intensify the high-frequency part of the jet noise spectrum. 


2. Scaling of aerodynamic noise at low Mach number 

An airflow of characteristic velocity U and length-scale L, with high enough 
Reynolds number pULIp (where p and p are the air’s density and viscosity), is a 
turbulent airflow. The chaotic sound field which it radiates through the 
surrounding atmosphere (with undisturbed sound speed c 0 ) is known as 
aerodynamic noise. 

This sound radiation — apart from any effects of solid boundaries (see 
below) — is precisely that which would be generated by quadrupole sources of 
strength T v per unit volume (Lighthill 1952, 1978), where T v stands for the 
difference between the momentum flux in the real airflow and that in a simple 
acoustic medium with sound speed c 0 . The most important term in T { \ is the 
convective flux pu t u, of a momentum component pu t carried by a velocity 
component Uj. 

Turbulent airflows at low Mach number U/c a are compact sources of 
aerodynamic noise because typical frequencies oi in the turbulence scale as V/L 
(Strouhal scaling) and therefore the compactness ratio <siL/c 0 is small (Lighthill 
1962). It means that differences in phases of emission for sounds reaching a 
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distant observer are small enough for the whole flow field to radiate effectively as 
a single source. 

This source may be of dipole type, with dipole strength F„ in cases when 
the vector F t represents a force acting between the turbulent airflow and a solid 
body immersed in it (Curie 1955, Ffowcs Williams & Hawkings 1969). Typically, 
F ( scales as pV*L 2 , so that its rate of change scales as pt/ 3 L; and then the radiated 
power 

<Fj I >/12jtpc* W 

scales as plPlJlc 2 '. a sixth-power dependence on flow speed. Also, the acoustic 
efficiency rj, defined as the ratio of radiated power to a rate of delivery 
(proportional to pU y L 2 ) of energy to the flow, scales as (t//c 0 ) 3 = A/ 3 . 

On the other hand turbulent airflows at low Mach number in the absence 
of any such solid body radiate effectively as a single quadrupole source, with total 
strength Q v scaling as plfL 3 (because the strength Ty per unit volume scales as 
pi/ 2 ). The acoustic intensity at a distant point, whose vector separation and scalar 
distance from the source are x, and r respectively, is then 

r ' 2 ) 2 >/ 1 6 k V p c 3 , (2) 

where the second time-derivative Qy scales as pU*L. Accordingly, the radiated 
power scales as plfOtc 0 S (an eighth-power dependence on flow speed; see 
Lighthill 1952, 1962) and the acoustic efficiency r\ scales as ( Ulc a ) 5 = M f 5 . 

At low Mach numbers, therefore, such aerodynamic noise radiation of 
quadrupole type is unimportant whenever dipole radiation due to fluctuating body 
force (with efficiency proportional to M 3 ) is also present (Crighton 1975). In other 
problems, however, such as the noise of a jet (with practically no fluctuating body 
force), the quadrupole source becomes dominant. For example, a turbulent jet of 
exit speed U radiates with an acoustic efficiency of order lO^CZ/c,,) 5 . 
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3. Aerodynamic noise at not so low Mach number 


The chaotic nature of turbulent flow implies that velocity fluctuations at points P 
and Q, although they are well correlated when P and Q are very close, become 
almost uncorrelated when P and Q are not close to one another. Here we recall 
that the correlation coefficient C between velocities u P and u Q is defined as 

C - <v, • v <? >/<vV>' i <rp>' 4 (3) 

in terms of the deviations v,, = u P - <u p > and v Q = u Q - <u Q > from their 
means. When two uncorrelated quantities are combined, their mean square 
deviations are added up: 

<v f * Vg> 2 ■ <v^> ♦ <v^> if C ■ 0 (4) 

(because the term 2<v P ’ v Q > vanishes). 

Theories of turbulence define a correlation length £, with C not far from 
1 (ii P and Uq well correlated) when PQ is substantially less than £ , and C not far 
from zero (u P and u Q almost uncorrelated) when PQ substantially exceeds t. 
Roughly speaking, different regions of size £ ("eddies”) generate uncorrelated 
sound fields, and the mean square radiated noise is the sum of the mean square 
outputs from all of the regions (Lighthill 1954, 1962). 

Typical frequencies cj in the turbulence are of order o> = v/£, where v is 
a typical root mean square velocity deviation <v*>*. For each region of size £ t 
therefore, the compactness condition that 0 )£/c o be small is satisfied when v/c a is 
small. 


Compactness, then, requires only that a root mean square velocity deviation 
v, rather than a characteristic mean velocity U, be small compared with c 0 . The 
associated restriction on Ulc a is less and can be satisfied at "not so low" Mach 
number. 
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On the other hand the sound radiated is no longer that of a single 
quadrupole source. Rather it is a combination of uncorrelated radiation patterns 
from different regions of size l , each with an intensity field (2) where 
while by equation (4) the intensity fields of different regions may simply be added. 
Therefore, on division by the volume t 5 of a region, we obtain the intensity pattern 
radiated per unit volume of turbulence as 

l 3 < (? v x < jr y r' J )*>/ 16 n J r J pc*. ® 


4. Doppler effects on frequency, volume and compactness 

Moreover the radiation from any single eddy is subject to a modification — the 
Doppler effect — as a result of the eddy being convected at "not so low" Mach 
number. A rather familiar element of Doppler effect on the radiation pattern of 
a moving source of sound is its shift in frequency, but there are changes also in its 
effective volume, and in its compactness. 

All of these Doppler effects on the sound received by an observer at a far- 
field location depend not on the speed V with which the source moves but on its 
velocity component >v in the direction of the observer. In the case of an observer 
located on a line making an angle 0 with a source’s direction of motion at speed 
V, this velocity component w has the form 

H- =* KCOS0 . (6) 

Then, while sound radiation of frequency travels a distance c a T during a 
single period T = 2ir/w, its source moves a distance wT nearer to the observer. 
Thus the wavelength X (distance between crests) is reduced to 
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k = c o T - wT = 27 i (c, -w)/« , 


(7) 


and the frequency heard by the observer (2?r divided by the time \/c 0 between 
arrival of crests)is increased to its Doppler-shifted value: the "relative" frequency 

= u(l - h ’(cj' 1 * o[l - Wcjcos&r 1 ; (8) 

though this, of course, may represent a decrease where the angle 9 is obtuse. 

The corresponding change in effective volume results from an effective 
change in the source’s dimension f in the direction of the observer (Lighthill 
1962). Because the near side N of the source region is closer by a distance l than 
its far side F, there is a certain "lag" r between the times of emission from F and 
N of sounds that reach the observer simultaneously. Then in the time t for sound 
from F to arrive it travels a distance cj , but the corresponding emission from N 
starts after a time lag r during which its distance from F in the direction of the 
observer has increased from t to ( + wr , and the sound emitted then travels a 
distance cjj - r). The condition for both sounds to reach the same point at time 
t may be written 

cj = t + wt + C 0 (t - x), giving t = tl(c 9 - w) . (9) 

The source’s effective dimension in the direction of emission is therefore altered 
to 

{ + vva = f(l = (uj co : (1®) 

a change by the same Doppler factor ujo) that modifies the frequency. Indeed, the 
eddy’s effective volume during emission is also increased by the Doppler factor 
w/o), because dimension in the direction of the observer is so increased whilst 
other dimensions are unaltered. 


Because the effective eddy volume f 3 appears to the first power in 
expression (5), whilst effective frequency occurs to the fourth power in the mean 
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square of the second time derivative f u , these Doppler effects produce an overall 
modification by five Doppler factors (Ffowcs Williams 1963, Lighthill 1963), 

(o,/w) 5 - (1 -w/cy - [1 -WOcosej* 5 , (11) 

in the intensity pattern (5) radiated per unit volume of turbulence. This gives a 
first indication of an important preference for forward emission from turbulence 
convected at "not so low" Mach number. 

On the other hand, it is essential to recognise how, as V/c 0 increases, 
Doppler effect tends also to degrade the compactness of aeroacoustic sources in 
relation to forward emission. Not only does a)f/c„ increase in proportion to Mach 
number, but an even greater value is taken by o) r (lc ot the ratio which must be small 
if convected sources are to be compact. A very marked restriction on the extent 
(1 1) of intensity enhancement for forward emission as Vlc a increases is placed by 
these tendencies (Ffowcs Williams 1963, Lighthill 1963, Dowling, Ffowcs 
Williams & Goldstein 1978). They can develop, indeed, to a point where the 
compact-source approximation (of low-frequency acoustics) may appropriately be 
replaced by its opposite (high-frequency) extreme: the ray-acoustics approximation. 
Thus, for supersonic source convection ( V/c„> 1), the relative frequency (8) 
becomes infinite in the Mach direction 

e = cos ; ( 12 ) 

moreover, it may be shown that radiation from the source proceeds along rays 
emitted in the Mach direction (Ffowcs Williams & Maidanik 1965). 

As this is a paper which devotes special attention to such radiation in the 
Mach direction, some immediate comments on its nature may perhaps be made. 
Equation (6) shows that the source’s velocity of approach w towards an observer 
positioned at an angle (12) to its direction of motion is the sound speed c a . On 
linear theory this means that different parts of a signal are observed simultaneously 
— which, of course, is the well known condition of stationary phase satisfied on 
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rays (Lighthill 1978). Sounds emitted by a source approaching at a speed w 
exceeding c a would be heard by the observer in reverse order (so that "pap pep pip 
pop pup" became "pup pop pip pep pap"!) but when w = c„ all the sounds (vowels 
and consonants!) would be heard together as one single "boom". 

From these preliminary comments the importance of two influences treated 
in later sections, that place limits on the signal propagated along rays, will already 
be clear. These are (i) the duration 5 of well-correlated emission from turbulent 
eddies and also (ii) nonlinear propagation involving a departure of the signal 
velocity from c a . 


5. Uniformly valid Doppler-effect approximations 

The correlation duration 5 for convected turbulent eddies is defined so that 
velocities at times separated by substantially more than 5 are almost uncorrelated 
while there is good correlation between velocities at times separated by 
substantially less than 5. This definition in terms of time separations (for the 
moving eddy) is directly parallel to the definition of i in terms of spatial 
separations. 

Combined use of correlation length ( and duration 5 affords an 
approximation to the radiation pattern from convected eddies that has some value 
at all Mach numbers. Thus it is a uniformly valid approximation, spanning the 
areas of applicability of the compact-source and ray-acoustics approximations. 

Figure 1 uses space-time diagrams where the space-coordinate (abscissa) 
represents distance in the direction of the observer. Diagram (a) for unconvected 
eddies approximates the region of good correlation, which must have spatial and 
temporal dimensions t and 5, as a simple smooth curve; actually, an ellipse with 
t and 5 as its axes. Diagram (b) shows how the region of good correlation is 
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sheared (sheared, in fact, by a distance w per unit time) for convected eddies 
approaching the observer at velocity w. 

Signals from a far point F and a near point N, in either case, reach the 
observer simultaneously — as do signals from other points on the line FN — 
provided that this line slopes by a distance c 0 per unit time. Then diagram (b) 
distinguishes three cases as follows: 

(i) in the compact-source case wlc 0 is small, and the space 

component of FN in diagram (b) is £(1 - w/c 0 )‘\ just as in 
equation (10) for the "usual" Doppler effect (which neglects 
the upper bound 5 on correlation duration); 

(ii) in the ray-acoustics case w/c 0 = I, and the space component of 

FN is c 0 5; 

(iii) in the intermediate case w/c a is but moderately less than 1, and 

the space component of FN is t multiplied by an 
enhancement factor 


[(1 - We.) 1 ♦ (t/c.Sfp • (13) 

Evidently, this form (13) of the enhancement factor comprehends all three cases 
and represents the effective augmentation of source volume due to convection. 

The enhancement factor (13) needs to be applied, not only to the volume 
term £ 3 in the quadrupole field (5), but also twice to each of the pair of twice- 
differentiated terms inside the mean square; essentially, because time- 
differentiations in quadrupole fields arise from differences in the time of emission 
by different parts of the quadrupole source region — and the time component of 
FN in diagram (b) is simply the space component divided by c 0 . As before, then, 
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five separate factors (13) enhance the intensity field; and, with w replaced by 
V cos 0, expression (11) for the overall intensity modification factor is replaced 
by 

{[1 - (K/c,)cosef + (t/c/) 1 }-* 2 . ( 14 ) 

This is a significant change wherever 1 - ( V/c a ) cos 0 is relatively small, and it 
tends to limit the predicted preference for forward emission (Ffowcs Williams 
1963, Lighthill 1963). 

It also gives an improved description of the influence of Doppler effect on 
the overall acoustic power output from convected turbulence. For example, the 
solid line in diagram (c) gives a log-log plot of the average value (spherical mean) 
of the modification factor (14) as a function of VI c a on the reasonable assumption 
that i = 0.6 V5. It will be noted that, as Vlc 0 increases, this average modification 
factor rises a little at first, but falls drastically like 5(V7c 0 )' 5 for V/c 0 significantly 
greater than 1 . 

On the other hand aerodynamic noise at low Mach number (see above) has 
an acoustic efficiency rj scaling as ( Ulc 0 ) 5 where U is a characteristic velocity in 
the flow. Actually, it would be permissible to take that characteristic velocity as 
the eddy convection velocity V ; although, if this were done in the case of a jet, it 
would be important to recognize that V is not the jet exit speed itself but takes 
values between 0.5 and 0.6 times the jet exit speed. Thus, an order of magnitude 
10 4 (C//c 0 ) 5 for i} in terms of jet exit velocity U corresponds to an order of 
magnitude 10' 3 (V'/c o ) 5 in terms of the eddy convection velocity V. 

The broken line in diagram (c) shows how this acoustic efficiency rj, 
supposed to take the value 10' 3 (V7c o ) J for small V/c 0 , is modified after 
multiplication by the average modification factor (solid line). This modification 
causes rj to approach a constant value of about 0.005 (aeroacoustic saturation) at 
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high Mach numbers, when slightly less than 1 % of jet power is radiated as 
aerodynamic noise. 

Supersonic jets are called "properly expanded" when they emerge — from 
appropriately shaped nozzles — as parallel flows. This is a contrasting case to that 
of a supersonic jet emerging as a non-parallel flow, which is necessarily followed 
by a sequence of shock waves (leading to augmentation of aerodynamic noise) in 
the characteristic "diamond" shock-cell pattern. 

Aeroacoustic saturation similar to that indicated in diagram (c) is observed 
for properly expanded supersonic jets, with acoustic efficiency becoming close to 
an asymptotically constant value of a little less than 10' 2 when the jet Mach number 
Ulc 0 exceeds about 2. These are the "extreme-speed" jets, with V itself (the eddy 
convection velocity) exceeding c„, for which aerodynamic noise is directed (Ffowcs 
Williams and Maidanik 1965) along rays inclined at the Mach angle (12). Before 
dedicating the remainder of the lecture to the nature of aerodynamic noise in this 
extreme-speed limit, I have attempted through the above discussion to exhibit that 
continuous trend of changing noise patterns which links it to the problem of 
aerodynamic noise at low Mach number. 


6. Nonlinear propagation of noise from extreme-speed jets 

Extreme-speed jets have just been defined as those properly expanded jets, with 
speeds more than twice the atmospheric sound speed c 0 , for which the eddy con- 
vection velocity V is itself supersonic and noise is primarily emitted in the Mach 
direction (12). Significant influences on this type of aerodynamic noise are 
exerted, not only by correlation duration (see above) for a supersonically moving 
eddy, but also by important consequences of nonlinear propagation. 
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Indeed, in a rather obvious analogy with a supersonically moving body, we 
may expect each eddy’s sound field not only (i) to be emitted in the Mach direction 
but also (ii) to take the form of a "supersonic boom" which, from nonlinear 
propagation effects, incorporates one or more shock waves. On the other hand, 
because it is a chaotic sequence of supersonically moving eddies that generates 
such waveforms, the jet’s near noise field must consist (as is, indeed, observed) 
of a thickly packed random assemblage of conical shock waves. 

Thus although, as just remarked, there may be some sort of continuous 
trend in the mechanisms of generation of aerodynamic noise between classical 
processes at lower Mach numbers and such extreme-speed jets, nevertheless studies 
of its propagation become quite different in this latter case. As a marked contrast 
to approaches that use the full nonlinear flow equations in the jet and simple linear 
equations of propagation outside it (generation being associated with differences 
between the equations), it becomes necessary to recognize important nonlinear 
effects on propagation itself. 

In the literature of aerodynamic noise for extreme-speed jets there already 
exist studies of the generation process that go far beyond the perfunctory sketches 
which I have included in this lecture. On the other hand, a properly detailed 
investigation of nonlinear effects on its propagation appears to have been neglected. 
For this reason, I have chosen to devote the rest of this lecture to such a detailed 
study of the nonlinear propagation of a thickly packed assemblage of conical shock 
waves once it has been formed. 

Briefly, this detailed study demonstrates how the inherent tendency of the 
shock strengths to decay (through conical spreading and internal dissipation) is 
counteracted in part by increases in strength that occur whenever adjacent shocks 
unite. A certain "bunching" tends to arise, because union of two adjacent shocks 
is found to increase the likelihood of further union with other neighbouring shocks. 
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The high-frequency part of the noise spectrum is made more intense than it would 
otherwise be by these bunching tendencies. 

A well known transformation of coordinates (Lighthill 1978) — suggested 
by the above-noted analogy with supersonic-boom theory (Whitham 1956) — is 
used (see Appendix) to reduce the problem of conical propagation, aimed at 
analysing how a temporal waveform varies with distance r along the Mach 
direction, to a plane-wave problem. In this latter problem the time t replaces r*, 
and attention is focussed on how a given spatial waveform varies with time. 

The spatial waveform assumed in this latter context is however immensely 
more complicated than the simple 'W-wave" form appearing in supersonic-boom 
theory. Instead, it is a spatially unlimited assemblage of random sawtooth waves 
of the type into which a general plane sound wave of large amplitude would evolve 
during a certain time t. 

This assemblage comprises randomly located shocks (with random 
strengths) separated by expansion waves in which the slope of the quantity "excess 
signal velocity" (excess, that is, over the undisturbed sound speed c 0 ) as a function 
of position takes the value 1/r. The inherent tendency of the shocks themselves is 
also to decay like 1 h\ but, as already described, this tendency is opposed whenever 
shocks unite with adjacent shocks. 

After an inverse transformation into the original variables describing conical 
propagation, t is replaced by r h but predicted shock strengths have to be divided 
by a further r' h factor. Thus the basic tendency of the shock strengths is now to 
decay like 1/r; but, yet again, the opposing tendency described above as 
"bunching" may prove to be important for the analysis of how aspects of 
aerodynamic noise are modified in the region of conical propagation that surrounds 
an extreme-speed jet. 
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7. Random sawtooth waves in transformed coordinates 


The classical transformation described in the Appendix addresses the problem of 
how a thickly packed random assemblage of conical shock waves will evolve with 
increasing r (distance in the direction of propagation) by reducing it to an 
interesting, yet hitherto neglected, problem of plane-wave propagation. This plane- 
wave problem, which will now be defined, treats the evolution in time of a 
spatially unlimited assemblage of random sawtooth waves. 

Because a random sawtooth wave is primarily important as a form of 
acoustic noise, any useful specification of such a wave must be one which 
facilitates identification of its noise spectrum. Now, noise spectra are determined 
in practice — whether in a physical or in a numerical experiment — from signal 
records of very great yet finite length which are Fourier analysed by F.F.T. 
techniques. 

Moreover any such Fourier analysis of a signal in an interval of great yet 
finite length expresses it essentially as a Fourier series, where successive terms in 
the series describe oscillations with very closely neighbouring frequencies. The 
noise spectrum is proportional to the squares of their amplitudes, regarded as a 
function of frequency, and the very close spacing between neighbouring 
frequencies allows it in practice to be depicted as a continuous curve. 

Actually, the Fourier series in question represents an exactly periodic 
function, with period equal to the length of the interval. Specifically, it represents 
the unique periodic function (with that period) which coincides with the signal 
record within the interval. 

These considerations lead us, in any study of spatially unlimited waves of 
random sawtooth type, to focus attention in practice on a random sawtooth wave 
which is specified on an interval of great but finite length L; while, for 
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completeness, its form outside that interval is defined by requiring it to be a 
periodic function with period L. There is no implication here that our interest is 
really confined to strictly periodic functions; on the contrary, we are interested in 
waves that, besides being spatially unlimited, have an everywhere random 
character. But we recognize that any physical or numerical experiment will 
confine attention to an interval of great but finite length L within which Fourier 
analysis of the signal gives a representation of it as an exactly periodic function of 
period L, and it is this representation which we find fruitful to analyse. 

Plane waves travelling in the ^-direction in a homogeneous medium with 
undisturbed sound speed c 0 may in nonlinear acoustics be conveniently described 
in terms of an independent variable 

X - x-c t t (15) 

and a dependent variable equal to the excess signal velocity 

v = u * c - c 0 , (15) 

where u is the fluid velocity and c the local sound speed (equal to c a + 0.2w in 
air). In a frame of reference moving at speed c ot as defined by the space 
coordinate X, any value of v is propagated at a velocity equal to v itself. 

This is the result which, whenever v is continuous, may be represented by 
the familiar partial differential equation 

dv/dt + vdv/dX = 0 , (I 7 ) 

with its characteristics of slope v in the ( t,X) plane. However, the system 
necessarily tends to develop shocks, which will here be treated as discontinuities. 
The speed of a shock is equal to the average of the smaller and greater values of 
v which appear ahead of and behind it; accordingly, the shock absorbs 
characteristics ahead of it by running into them, while absorbing characteristics 
from behind as they mn into it. 
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It may be shown (for example, by differentiating (17) with respect to X) that 
the reciprocal slope ( dv/dX )'* of continuous parts of the wave must increase at unit 
rate along any characteristic. It means that any negative value of the reciprocal 
slope must after a finite time increase to zero, corresponding to infinite slope. 
This, of course, is when a shock appears, with its subsequent propagation governed 
by laws quoted above. 

By contrast, any positive initial value of the reciprocal slope must grow 
indefinitely at unit rate, which after a time t adds a term t to that initial value. 
When t has become large, this added term t has become dominant over the initial 
value (at least, if this is not too big), so that to a close approximation the slope 
dv/dX itself takes the value 1/r. Indeed, a classically familiar argument why this 
must prove to be rather a good approximation is that, when t has become large, 
all characteristics with relatively bigger initial values of the reciprocal slope have 
disappeared through running into shocks; essentially, because of their tendency to 
be close to characteristics with negative initial values of reciprocal slope (from 
which shocks necessarily develop). The continuous parts of the waveform are 
dominated, therefore, by characteristics on which the slope has become close to 
1/r. 


These are the reasons why an initial random acoustic wave of large 
amplitude develops after time r into what we define as a random sawtooth wave. 
This comprises shocks at positions 

X = x n (- « < „ < oo) (18) 

with X n a decreasing function of the integer n\ where the values of v behind and 
ahead of the nth shock are 

v B and V„ - , (19) 

and where we shall describe as the "strength" of the mh shock the discontinuity z n 
between the values of v behind and ahead of it. Moreover, the distribution of v 
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has the slope 1 It between shocks, so that the value v„ - z H of v just ahead of the nth 
shock is related to its value v „ . , just behind the (n - l)th shock by the equation 

v, -V, - X n yt . (20) 

The decision to assume the wave periodic with period L, which has already 
been carefully explained, has certain consequences for the quantities (18) and (19). 
At any particular time t there will be a certain specific number N of shocks within 
a single period; here, N is a function of the time, taking integer values, which is 
reduced discontinuously by 1 at every instant when a union of two shocks occurs. 
The periodicity then implies that 

\. N “ v „ • Vw = > x *.n - x . " 1 • W 

so that equation (20), summed between n = 1 and n — N, gives 

E *. = ■ < 22 > 

n ■ I 

This equation, indicating a balance between the net compressive and expansive 
effects in the sawtooth wave, ensures that the period L remains unchanged because 
the sum of the shock strengths on the left-hand side is found (see below) to vary 
as lit. 


Figure 2 illustrates the form of a random sawtooth wave at a particular time 
t = T, with the number N of shocks inside the period 0 < X < L equal to 25. 
The randomness of this wave derives from the fact that a random number generator 
was used to determine all the strengths z n of the 25 shocks and all the spacings 
~ %n-i - Xn between shocks, subject to 

(i) the need, which equation (21) demonstrates, for the sum of the spacings 

h n to be L; 

(ii) the corresponding condition (22) on the sum of the strengths; and 

(iii) restrictions on each value of h n or of tz„ to lie between 0.01L and 

0.09L. 
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Thus the spacings and strengths, subject to (i) and (ii), take random values between 
these upper and lower bounds. (Here, while the upper bound is intended to reflect 
some sort of limitation on sizes and strengths of noise sources, the lower bound is 
imposed mainly for convenience of graphical representation.) The above values 
fix uniquely a plot of v against X with zero integral over a period, and the 
evolution of this particular sawtooth wave is exhibited in section 11 below. 


For a random function v of X defined in an interval 0 < X < L (here, a 
particular single period for the function), its noise spectrum as defined so that 
P(k)dk represents the contribution to the function’s mean square from wavenumbers 
between k and k + dk is classically given as 

P(k) = < jj^ve'^dxj 2 > (23) 


where the angle brackets denote a statistical mean value. The high-wavenumber 
behaviour of this noise spectrum, for a random sawtooth function v dominated by 
discontinuous changes z n at points X n as in (18) and (19), is 
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(24) 


so that it depends on the sum of the squares of the shock strengths z n . 


It is this dependence which makes unions of shocks important. At any such 
union, the sum of the shock strengths themselves remains unchanged; however, the 
sum of their squares is necessarily increased, so that the high-frequency part of the 
noise spectrum is intensified. 


8. Sawtooth-wave evolution up to when a first union occurs 

The dynamics of the shocks is analysed next up to when a union first occurs. At 
each instant the basic law 
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*. - V. - vfe. . 


(25) 


which expresses the velocity X n of a shock wave as the average of the values (19) 
of v behind and ahead of it, implies also that 

because, during a small time the shock absorbs into itself from behind a section 
'/iz n bt of smooth waveform with slope \!t. 


Now careful study of the system of equations (25), (26) and (20) shows that 
their completely general solution takes the form 


Y '. \ + 0720 , X. 


+ 'AY. * bj 


(27) 


where a n , b n and Y„ are constants satisfying 

v, » 8 > 

This solution gives a constant value Y n to the product of the shock strength z„ with 
the time t. Here, as noted earlier, t is measured from an origin (r = 0) when the 
waveform was initiated, so that it later developed into a sawtooth wave with the 
slope 1/r for all smooth sections thereof. 


Great interest attaches to the distance h n = X n _, - X n between two adjacent 
shocks. Equations (27) and (28) show that 

K - + r ») + (*-» " b *Y ■ (29) 

Now, as noted earlier, the analysis in this lecture is only concerned with 
propagation of random sawtooth waves once they have appeared; say, from t = T 
when h n takes the form 
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*■ = v*(r,_| ♦ r,) jgr. 


(30) 


Eliminating the unknown constant b n .j - b n between (29) and (30) we obtain an 
equation for the evolution of h n in the form 

K - «*>[*. - (n-i + *.)<>} . (31) 

where it proves useful to introduce a new time-like variable a. The equations 

a • (t - Dpt , t - 7/(1 - 2o) (32) 

make a an increasing function of t and map the unbounded interval T<t<eo into 
the finite interval 0 < a < 'A. It turns out that a is a highly convenient measure 
of the times at which unions of shocks occur. 


As just a preliminary example of this, we note from equation (31) that the 
very first union of two shocks occurs when 


a 




(33) 


since it is this value of a that first allows one of the h n (distances between shocks) 
to fall to zero. Because of the periodicity assumption which was exhaustively 
discussed earlier, the minimum (33) must of course be attained not only for a 
particular value n within the period 1 < n ^ N but also for the corresponding 
value ( n ± mN, with m an integer) in any other period. Attention is here 
focussed, however, on a single period — except for the fact that when equation 
(28) or (33), or any other equation involving n - 1, is applied to (say) n = 1, 
periodicity is used to interpret n - 1 as N, 


On the other hand, there is no need in a random sawtooth wave to consider 
the possibility (really an impossibility — since it would occur with zero measure 
in a probability space of random variables) that the same minimum value (33) 
might be attained for two different n in the period 1 <, n <, N. It will be assumed 
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rather that the first union occurs for a particular n and a, as specified by equation 
(33). Also, because equations (22) and (27) imply that 

N N KA\ 

E >■, ■ i - E «. (34) 

the minimum (33) is readily shown to satisfy the requirement, a < 'A, for it to 
correspond to an actual time t > T for first union between shocks. 

Until this time, of course, every Y n (defined as the product of the strength 
of the mh shock with the time t ) has remained constant. Then a first union — say, 
of the mh and ( n - l)th shocks — produces addition of their strengths z n and z„. h 
with consequences noted earlier for the high-frequency form (24) of the noise 
spectrum. 

In the meantime, there is a reduction by 1 in the number N of shocks in a 
single period, but the sum on the left-hand side of (34) continues to take the value 
L simply because the new shock, replacing those associated with the T-values Y n 
and Y m .„ takes on their added T-value Y n + Y n _,. Similarly equation (34) remains 
unchanged at all later unions — other properties of which will now be investigated. 


9. Formulas specifying all later unions of shocks 

In this section, the tendency to "bunching" of shocks — that is, an increased 
likelihood of union between shocks that have already participated in union with 
other shocks — is quantified by means of a general formula for the time of union 
of two shocks which allows for all preceding unions. The two shocks to be 
considered are taken as (i) a shock originally numbered n, into which other shocks 
have merged from behind, and (ii) a shock originally number n - 1, which has run 
into (and merged with) various shocks ahead of it. 
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Evidently, this is the most general possible merger. Indeed, when any pair 
of shocks unite, it is permissible to "identify" the back one with an original shock 
(numbered n, say) into which others have merged from behind, and the front one 
with an original shock (numbered n - 1 , say) that has run into others ahead of it. 

Specifically, we shall suppose that the shock originally numbered n has 
participated in a total number m n of unions with shocks from behind, with T-values 

Y Bm , occurring at times t ■ ■ 1 to m,) . (35) 

Similarly, we assume that the shock originally numbered n - 1 has participated in 
a total number m„., of unions with shocks ahead of it, with T-values 

Y AM , occurring at times t * t Am [m 9 1 m H .\) ■ (36) 

It is convenient to embark upon this problem by considering first the 
dynamics of the shock originally numbered n immediately after it has undergone 
just the first of the above unions, at t = t Bl . Because we "identify" this merged 
shock with the shock originally numbered n, we continue to use the subscript n in 
quantities such as (18) and (19) associated with this shock. Then the differential 
equations (25) and (26) continue to apply to this united shock, as does the 
relationship (20), so that the form (27) of this system’s general solution is 
unaltered; However, the constants a n , b n and Y„ are changed to new values 

a. ■ y„ . K • , r. • r„ ; (37) 

while (28) is still satisfied since a n . x has not been changed. 

Out of these new values (37) for a n , b H and Y n , the last (added K-values for 
uniting shocks) has been explained above. Also, we may readily verify that the 
new values for a n and b n are the only ones which satisfy two essential conditions: 
that union at time t = t BX makes no discontinuous change in the position X n of the 
nth shock, while increasing v„ (the value of v behind it) by z Bl = Y Bl /t BX (the 
merging shock’s jump in v). 
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Similar studies show that, in each subsequent union of the nth shock with 
shocks from behind as specified in (35), the constants receive further changes as 
in (37) but with and t ^ replacing Y BX and t BX \ so that, after all of them, a n , b„ 
and Y n have become 

». - E I'm . K * E • r . *£>■..' (38) 

M-l *»■! 

Other careful studies of the unions of the ( n - l)th shock with shocks ahead of it 
as specified in (36) demonstrate that, when all are completed, a n _ ,, b n . x and Y„. x 
have become 


'■-I 


E !■(. • V, - E (Wt.) ■ r.-x * E K 


Am 


(39) 


These results allow both X n and X n _ x to be obtained by use of (27), so that 
the separation h n = X n . x - X n between the shocks can be expressed as 


*. = + y «) + - b -Y - 1/2 E Wa* - 1) - Yb.V'b- - l ) • 


(40) 


HI") 


91*1 


Here, equation ( 30 ) may be used to substitute for b nA - b n in terms of the value H n 
of the separation h n at the initial time t = T. Then, with times t, t ^ and t ^ 
substituted in terms of the corresponding cr-values from (32), we obtain 

(1 - 2 o)h„ = //, - (Y^ ♦ Y n )o - "f YjJa - a,.) - £ Y^a - o ta ) . (41) 

« * 1 M-i 

Equation (41) gives an elegant and completely general expression for the 
value of the time-like variable a when the nth and (n - l)th shocks unite (so that 
h n becomes zero). The form of this expression becomes particularly instructive if 
we use a n to signify the ratio 

= Wn-l + Y ') • (42) 


Equation (31) shows how, provided that o a < Vi, this value of a specifies the time 
at which the nth and (n - l)th shocks would have united if no other unions of 
shocks had occurred first. 
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On the other hand, the actual value of a corresponding to the time at which 
the shocks unite is given, by using the substitution (42) after putting h n = 0 in 
(41), as 


o 


PVi - ♦ E r*** + E 


JH-I 


M-l 


♦ r, * E y a. * E >» 

M-l 


(43) 


Expression (43) may be recognized as a weighted mean of the a-values <x„, a ^ and 
cr^ with weights given in terms of shock strengths according to the following 
simple rules. 


The weight attached to a Am , or to o^, is K^or respectively; in other 
words, it is the y-value (strength multiplied by t) of the shock that merges when 
a = with the (n - l)th shock or when a = a Bm with the mh shock. By 
contrast, the weight attached to a n is the initial combined y-values of the mh and 
(n - l)th shocks which would have united at o - o n (provided that cr n < 'A, though 
otherwise not at all) if no preceding unions had occurred. The sum of all these 
weights, which constitutes the denominator of (43), represents the overall y-value 
of the united shock that results from union of the mh and ( n - l)th shocks after 
they have respectively undergone all the mergers (35) and (36). 


10. Bunching tendencies and their implications 

Very simple properties of the weighted mean (43) tell us that, because all these 
mergers at the times (35) and (36) have preceded the union of the mh and ( n - l)th 
shocks, so that 

o > c Am and o > a ^ for all m , (44) 

therefore the quantity a„, which is the only other among all the a-values of which 
(43) represents the weighted mean, must satisfy the opposite inequality a < o n . 
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In words, the union of the nth and ( n - l)th shocks takes place earlier, in 
consequence of any preceding unions in which either have participated, than would 
have been the case if no such preceding unions had occurred. 

But still more valuable than any such purely qualitative statement is the 
quantitative expression (43) for the reduced value of a at the instant of union. For 
example, in cases when the Y-values of the nth and ( n - l)th shocks were initially 
rather small — before they merged with stronger shocks — it is of course quite 
possible for a n , as defined by the ratio (42), to be considerably greater than 'A. 
This implies that the shocks are initially moving farther apart, so that their union 
could never have taken place if neither had merged with other shocks. Yet, in the 
weighted mean (43), the combined weights attached to all the cr^ and cr^ may 
greatly exceed that attached to a n if the Y - values of the merging shocks amount in 
total to much more than the initial combined value + Y n \ and, evidently, this 
allows the value (43) for a (associated with the time of union of the mh and 
(n - l)th shocks) to be only a little bigger than the greatest of the cr Am and <r Bm . 

Such considerations allow us to infer that a random sawtooth wave may be 
subject to a very marked "bunching" process. This is a process of "snowball" type 
where early local unions of shocks act to stimulate further unions with 
neighbouring shocks. Because this process, while leaving the sum of all the 
T-values unchanged, produces an increase in the sum of their squares, it may have 
an important effect of enhancing (as expression (24) shows) the high-frequency part 
of the noise spectrum. 


II. An algorithm facilitating numerical experiments on bunching 

Another useful route to the quantitative study of the bunching process is through 
numerical experiments. In principle, these might be attempted by computing 
"weak" solutions of equation (17) and scrutinizing their evolution from instant to 
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instant to identity those times when unions of shocks occur. On the other hand 
such identifications can be a little inconvenient and somewhat lacking in precision 
with practical values of the grid spacing. 

A much more straightforward approach is one that treats all shocks as 
discontinuities and which makes use of just a single time-step between each union 
and the next. The basic result (33), for defining when an initial set of N shocks 
within a period L becomes reduced to a set of N - 1 shocks after the first union of 
two shocks has occurred, may be applied successively for this purpose in the 
following algorithmic treatment. 

The algorithm starts from N shocks filling the period L at a given time t, 
with T-values Y n (n = 1 to AO and with spacings h n = X n ., - X n between shocks, 
where 

£ h„ - L and, by (22), £ Y K - L . (45) 

N«l "■! 

The form of the algorithm is a respecification of all these variables at the instant 
when the first union of any pair of adjacent shocks has occurred. (This algorithm 
can then, systematically, be given successive applications to determine subsequent 
developments, including bunching.) 

The most obvious change in one of the variables after that union is that N 
has been replaced by W - 1. All of the other changes, however, depend critically 
on those quantities which have been shown to define when a union first occurs; 
namely the minimum (33) and the value of n for which the minimum is attained. 

In the present case we define, then, 

«„ ■ [W-i * >'.)] • «* 1 “ d ** • - " • (4<) 

and note that, according to equations (45), this minimum value cannot exceed 'A. 
Then the time t needs to be respecified as the time 
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</(l - 2<g 


( 47 ) 


at which a union first occurs. At that instant the /nth and (m - l)th shocks have 
united, and are designated thereafter as the (m - l)th, with y-value Y mA 4- Y m . For 
other values of n, the nth shock retains its T-value Y„, and is still numbered as the 
nth shock for n < m - 1 but now needs to be renumbered as the (n -l)th shock for 
n > m. Also, equation (31) tells us that the spacing h n is changed to 

(*. - * *•.)]/(> - 2 «.) <48) 

for n < m - 1. The same change occurs for n > m + 1 except that the value of 
n is simultaneously replaced by n - 1 . 

A special note is required to the effect that quantities involving Y nA need 
when n - 1 to be reinterpreted (by periodicity) with Y nA replaced by Y N . Similarly, 
when m = 1, the united shock is then renumbered as the (N - l)th shock, with 
7-value y n + y,. 

The overall algorithm may be summarised as follows. Given a positive 
integer N and a time t, and a set of intervals and of T-values (shock strengths 
multiplied by the time) h„ and Y n for n = 1 to N, satisfying equations (45), we 
determine the minimum (46) where Y„ is to be interpreted as Y N . Then the 


algorithm 

replaces each of 

by the new value 

(a) 

N 

N - 1 

(b) 

t 

m - 2o 

(c) 

Y n for n < m - 1 

Y n 

(d) 

Y m -i if'" > 1 

Y m . i + Y m 

(e) 

Y na if m = 1 

y n + y, 
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(f) 

Y„ for m ^ n ^ N - 1 
(unless n = N - 1 and m = 1) 

n.i 

(g) 

h„ for n < m 

ik - o-.cn.. + n>i/< i - 2 o.) 

(h) 

h n for m £ n ^ N - 1 

[*»*. - o.cn + n«Mi - 2oj 


but it should perhaps be mentioned that, evidently, there are no values of n for 
which (c) can be applied if m - 1 or 2, and none for which (0 can be applied if 
m = N. Similarly, there are no values of n for which (g), or else (h), can be 
applied if m = 1 , or m = N, respectively. Subject to these interpretations, the 
algorithm can be very readily executed, and then re-executed as many times as 
required to determine all subsequent developments of the random sawtooth wave. 

Figure 3 shows the results of applying this algorithm to the random 
sawtooth wave illustrated in Figure 2. That configuration at the initial time t = 7, 
with its twenty-five shocks, occupies the bottom of the diagram; the rest of which 
depicts the paths of the shocks, including all unions, for T < t <, IT. 

During this period there have occurred one bunching of nine shocks, two 
bunchings of five shocks and one bunching of three shocks; on the other hand, 
three out of the original twenty-five shocks have avoided participating in any 
unions. All of this appears consistent with the suggested definition of bunching as 
a tendency "for shock waves that have already formed unions with other shock 
waves to acquire an increased proneness to form further unions". 

It may be natural to ask whether Figure 3 represents in any sense a 
"selected" example; this question, however, has a negative answer. Figure 2 was 
produced, as explained in section 7, by means of a random number generator, and 
the algorithm of this section was simply applied to the very first wave so 
generated. I shall be quite content to leave to later investigators the pleasure of 
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executing all those numerous "runs" of the algorithm that may be necessary to 
establish statistical trends! 

At each and every union, as already mentioned, the sum of the Y n remains 
constant but the sum of their squares increases discontinuously. Figure 4 shows this 
process in the case of the sawtooth waveform evolution illustrated in Figure 3. 
Within the interval T < t < IT quite a pronounced rise in E^, 2 is found, and it is 
intriguing that this discontinuously increasing dependence on t deviates only a little 
from a straight line through the origin, which would represent a simple 
proportionality to t. 

The implication here is that, although the inherent tendency of each shock 
strength z n = A/' is to vary as t A , nevertheless the sum E z„ 7 which appears in the 
high-frequency noise spectrum (24) shows a variation very different from r 2 and, 
indeed, one rather close to The bunching process, in short, significantly 
enhances the high-frequency part of the noise spectrum. 

12. Conclusion 

This has been a Lecture which, in the words of its title, has aimed at exploring just 
"some aspects" of jet noise fields. Indeed, even that single aspect which has been 
covered in most detail — the nonlinear propagation of the noise signal from 
extreme-speed jets — has been investigated only in a noise field quite close to the 
jet, where a rather consistently conical propagation of noise is to be found. This 
noise has its origin in the jet mixing region, separating an internal core flow of 
velocity U from the outside atmosphere, with eddies in this region convected at an 
approximately constant velocity V which defines the Mach direction (12) for noise 
emission. 

On the other hand, there exists only a finite length of jet mixing region, 
beyond which the core flow vanishes and eddy convection velocities immediately 
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start to fall off with distance. The beam of truly conical noise propagation, then, 
is only of finite width and, accordingly, can maintain its conical character for only 
a limited distance from the jet axis. Beyond that distance the propagation must 
change progressively to a propagation which — although still directional — is 
beginning to exhibit an essentially spherical attenuation. 

Nonetheless, although conscious of all this, I have here placed emphasis on 
the field of conical noise propagation rather close to the jet for two main reasons. 
The first applies if our concern really is with just an isolated jet radiating into still 
air. 


Then general nonlinear acoustics (Lighthill 1978) tells us that waveform 
modifications due to nonlinear propagation effects are enormously greater in a 
conically spreading wave than in a spherically spreading wave; specifically, the 
square-root transformation (A4) derived in the Appendix below is replaced by a 
logarithmic transformation for a spherical wave which, accordingly, undergoes far 
less modification from nonlinear propagation. The conical field is specially 
important, then, as the region of occurrence of most of the spectrum-modification 
effects of nonlinear propagation. 

But theoretical studies may have even more importance for the problem of 
how aero-engine installations are to be designed so as to reduce ground levels. In 
this context a jet noise field where propagation is directionally concentrated (say, 
in the Mach direction) may offer engine-installation designers a rather special 
opportunity to incorporate shielding devices aimed at limiting community noise. 
This, then, is yet another good reason why, in this inaugural Theodorsen Lecture, 
I have presented to my distinguished audience quite a detailed analysis of how the 
noise spectrum in that conical field with its "random assemblage of shocks" may 
be influenced by the general tendency on which I have affectionately ventured to 
bestow the name "bunching". 
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Appendix. Conical-wave problem transformed into plane-wave case 

A classical transformation converts the problem of conically propagating waves 
incorporating relatively weak shocks into a plane-wave problem. The 
transformation is effective in any region which for linear acoustics would be 
described as a far field, with amplitudes varying as r* along linearised 
characteristics. 

The corresponding law for nonlinear acoustics, in continuous parts of the 
field, is that amplitudes vary as r' A along characteristics as defined for the 
nonlinear problem (Whitham 1956). In both cases r represents propagation 
distance measured from a zero on the axis of the cones, but a certain lower limit 
for r exists beyond which this far-field representation can be applied. 

In the conical-wave problem treated here, it may be convenient to use a 
notation t c to signify the time, to avoid confusion with the use of t as the time for 
the plane-wave problem, treated above, into which the conical-propagation problem 
will be transformed. Moreover a suitable dependent variable (also designated by 
subscript c) will be the defect in reciprocal signal velocity (Lighthill 1978), 

v £ “ c p ‘ - (“ + c)' 1 . (A*) 

rather than the excess (16) in the signal velocity itself. 

Now the condition that this dependent variable varies as r' A along 
characteristics satisfying 

dt c ■* (u + c)' ] dr = (c;‘ - v^jdr (A2) 

may be written 

[djdr + (c; 1 - v c )dldt c )(r\) = 0 . (A3) 

Then a simple transformation of coordinates 

2 r* - t , c~ l r - t c = X ( A4 ) 

transforms the operator in braces into 
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(A5) 


r^d/at * c' 0 'didx - (c;‘ - vja/ax ; 

so that, if the dependent variable is changed also to v = r\, the equation for v 
becomes 

(d/dt * vd/dX)v * 0 , (A6) 

exactly as in equation (17) above. 

It follows that all the results for the plane-wave case can be directly applied 
to the conical-wave problem, with v standing for r A v c and t for 2 r h . The temporal 
evolution (as t increases) of the spatially unlimited waveform (in <. X < 00 ) 
for plane waves is then seen, according to (A4), as a representation of the spatial 
evolution (as r increases) of the temporally unlimited waveform (in -oo < t c < <») 
for conical waves. 

In this representation the important quantity Y n , which remains constant for 
each shock until it unites with another shock (when their K-values add), stands for 
the product of / with the change Av between the values of v ahead of and behind 
the shock. Its corresponding meaning for conical propagation is therefore 2rAv f , 
where Av c is the shock strength z„. Accordingly, in a conically propagating 
random sawtooth wave, the inherent tendency of shock strengths (while their 
T-values remain constant) is to decay like the reciprocal r' 1 of the propagation 
distance r from the axis of the cones, although this tendency is opposed wherever 
shocks unite (through their K-values adding). 

In the meantime the second of equations (A4) implies that the quantity h n , 
which in the plane-wave case stands for the distance between adjacent shocks at a 
given time, represents in the conical-wave problem the interval between times of 
passage of adjacent shocks past a given position with a specified value of r. Also 
the quantity T, used above to signify the time when a plane wave of random 
sawtooth type can be considered established, must in the conical-wave problem be 
given the value 2R' A , where R stands for the corresponding propagation distance 
from the axis at which such a wave has become established. 
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Calculations of times for unions of shocks in the plane-wave case treated 
above are expressed in terms of a nondimensional variable a, whose defining 
equations (32) correspond for conical propagation to the equations 

a - (r* - **)/ 2r* , r - R/( 1 -2 o) 2 . ( A7 ) 

Expressions for the values of a at which shocks unite, like (33) for the first such 
union or (43) for a quite general subsequent union — as given in terms of the 
T-values of shocks and the initial values H n of h n (now the initial time intervals 
between shocks) — can then immediately be interpreted, by (AT), as specifying the 
distance r required for merger of adjacent shocks. This in turn can be used to 
quantify tendencies to bunching in conical-wave propagation. 

Finally, if numerical experiments like that of section 1 1 indicate a general 
tendency for LY* to increase in approximate proportionality to / = 2 r' A , this has 
implications for the quantity 

E *.' ■ E <rjv? < A8 > 

*•1 fl“l 

which equation (24) associates with the high-frequency part of the noise spectrum. 
Although the inherent tendency of z n to decay at r 1 would cause this sum (A8) to 
vary as r 2 , bunching tendencies modify this decay law for high-frequency noise 
into a variation as r in . 
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Figure Captions 

Figure 1 . A uniformly valid Doppler-effect approximation. 

Diagram (a). Space-time diagram for unconvected "eddies" of correlation 
length i and duration 8. 

Diagram (b). Case of "eddies" convected towards observer at velocity 
w = V cos 0, being Diagram (a) sheared by a distance w per unit 
time. Here, straight lines sloping by a distance c 0 per unit time 
represent emissions received simultaneously by the observer. 

Case (i): w/c 0 small. Case (ii): w/c 9 = 1. Case (iii): intermediate 
value of w/c 0 . 

Diagram (c). Solid line: average (spherical mean) of the modification 
factor (13). Broken line: acoustic efficiency, tj, obtained by applying 
this factor to a low-Mach-number "quadrupole" efficiency of (say) 
10- 3 (V/c„) 5 

Application to jets. For a jet of exit speed (/, a typical eddy convection 
velocity V takes values between 0.5U and 0.6 U. In order of 
magnitude terms, then, the efficiency jj makes a transition, at about 
U/c 0 = 2, between values around 10 ^(U/c„y and an asymptotically 
constant value of a little less than 10' 2 . 

Figure 2. A random sawtooth wave, with the number N of shocks inside the 

I 

period 0 < X < L equal to 25; and with their spacings and strengths 
determined, subject to conditions (i), (ii) and (iii), by a random number 
generator. 

Figure 3. Results of applying the algorithm of section 11 to the initial waveform 
of Figure 2. This is reproduced at the bottom of the diagram; the rest of 
which depicts (by plots of XIL against t/T) the paths of the shocks, 
including all unions, for T < t £ IT. Reading from left to right, note two 
bunchings of five shocks each, one of three shocks and one of nine shocks, 
interspersed with three shocks that have avoided participating in any unions. 
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Figure 4. A plot of £ (Y a /L)* against tIT, for the sawtooth waveform evolution 

Jl-1 

illustrated in Figure 3, is used to indicate how bunching enhances the high- 
frequency part of the noise spectrum. 
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